Progenitor identification and SARS-CoV-2 infection in human distal lung organoids.

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange. Three-dimensional in vitro human distal lung culture systems would strongly facilitate the investigation of pathologies such as interstitial lung disease, cancer and coronavirus disease 2019 (COVID-19) pneumonia caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here we describe the development of a long-term feeder-free, chemically defined culture system for distal lung progenitors as organoids derived from single adult human alveolar epithelial type II (AT2) or KRT5+ basal cells. AT2 organoids were able to differentiate into AT1 cells, and basal cell organoids developed lumens lined with differentiated club and ciliated cells. Single-cell analysis of KRT5+ cells in basal organoids revealed a distinct population of ITGA6+ITGB4+ mitotic cells, whose offspring further segregated into a TNFRSF12Ahi subfraction that comprised about ten per cent of KRT5+ basal cells. This subpopulation formed clusters within terminal bronchioles and exhibited enriched clonogenic organoid growth activity. We created distal lung organoids with apical-out polarity to present ACE2 on the exposed external surface, facilitating infection of AT2 and basal cultures with SARS-CoV-2 and identifying club cells as a target population. This long-term, feeder-free culture of human distal lung organoids, coupled with single-cell analysis, identifies functional heterogeneity among basal cells and establishes a facile in vitro organoid model of human distal lung infections, including COVID-19-associated pneumonia.

Progenitor identification and SARS-CoV-2 infection in long-term human distal lung organoid cultures.

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange and is affected by disorders including interstitial lung disease, cancer, and SARS-CoV-2-associated COVID-19 pneumonia. Investigations of these localized pathologies have been hindered by a lack of 3D in vitro human distal lung culture systems. Further, human distal lung stem cell identification has been impaired by quiescence, anatomic divergence from mouse and lack of lineage tracing and clonogenic culture. Here, we developed robust feeder-free, chemically-defined culture of distal human lung progenitors as organoids derived clonally from single adult human alveolar epithelial type II (AT2) or KRT5 + basal cells. AT2 organoids exhibited AT1 transdifferentiation potential, while basal cell organoids progressively developed lumens lined by differentiated club and ciliated cells. Organoids consisting solely of club cells were not observed. Upon single cell RNA-sequencing (scRNA-seq), alveolar organoids were composed of proliferative AT2 cells; however, basal organoid KRT5 + cells contained a distinct ITGA6 + ITGB4 + mitotic population whose proliferation segregated to a TNFRSF12A hi subfraction. Clonogenic organoid growth was markedly enriched within the TNFRSF12A hi subset of FACS-purified ITGA6 + ITGB4 + basal cells from human lung or derivative organoids. In vivo, TNFRSF12A + cells comprised ~10% of KRT5 + basal cells and resided in clusters within terminal bronchioles. To model COVID-19 distal lung disease, we everted the polarity of basal and alveolar organoids to rapidly relocate differentiated club and ciliated cells from the organoid lumen to the exterior surface, thus displaying the SARS-CoV-2 receptor ACE2 on the outwardly-facing apical aspect. Accordingly, basal and AT2 apical-out organoids were infected by SARS-CoV-2, identifying club cells as a novel target population. This long-term, feeder-free organoid culture of human distal lung alveolar and basal stem cells, coupled with single cell analysis, identifies unsuspected basal cell functional heterogeneity and exemplifies progenitor identification within a slowly proliferating human tissue. Further, our studies establish a facile in vitro organoid model for human distal lung infectious diseases including COVID-19-associated pneumonia.

Advances in Proximity Ligation in situ Hybridization (PLISH).

Understanding tissues in the context of development, maintenance and disease requires determining the molecular profiles of individual cells within their native in vivo spatial context. We developed a Proximity Ligation in situ Hybridization technology (PLISH) that enables quantitative measurement of single cell gene expression in intact tissues, which we have now updated. By recording spatial information for every profiled cell, PLISH enables retrospective mapping of distinct cell classes and inference of their in vivo interactions. PLISH has high sensitivity, specificity and signal to noise ratio. It is also rapid, scalable, and does not require expertise in molecular biology so it can be easily adopted by basic and clinical researchers.

Single-cell Wnt signaling niches maintain stemness of alveolar type 2 cells.

Alveoli, the lung's respiratory units, are tiny sacs where oxygen enters the bloodstream. They are lined by flat alveolar type 1 (AT1) cells, which mediate gas exchange, and AT2 cells, which secrete surfactant. Rare AT2s also function as alveolar stem cells. We show that AT2 lung stem cells display active Wnt signaling, and many of them are near single, Wnt-expressing fibroblasts. Blocking Wnt secretion depletes these stem cells. Daughter cells leaving the Wnt niche transdifferentiate into AT1s: Maintaining Wnt signaling prevents transdifferentiation, whereas abrogating Wnt signaling promotes it. Injury induces AT2 autocrine Wnts, recruiting "bulk" AT2s as progenitors. Thus, individual AT2 stem cells reside in single-cell fibroblast niches providing juxtacrine Wnts that maintain them, whereas injury induces autocrine Wnts that transiently expand the progenitor pool. This simple niche maintains the gas exchange surface and is coopted in cancer.

Automated cell-type classification in intact tissues by single-cell molecular profiling.

A major challenge in biology is identifying distinct cell classes and mapping their interactions in vivo. Tissue-dissociative technologies enable deep single cell molecular profiling but do not provide spatial information. We developed a proximity ligation in situ hybridization technology (PLISH) with exceptional signal strength, specificity, and sensitivity in tissue. Multiplexed data sets can be acquired using barcoded probes and rapid label-image-erase cycles, with automated calculation of single cell profiles, enabling clustering and anatomical re-mapping of cells. We apply PLISH to expression profile ~2900 cells in intact mouse lung, which identifies and localizes known cell types, including rare ones. Unsupervised classification of the cells indicates differential expression of 'housekeeping' genes between cell types, and re-mapping of two sub-classes of Club cells highlights their segregated spatial domains in terminal airways. By enabling single cell profiling of various RNA species in situ, PLISH can impact many areas of basic and medical research.

Ultrasensitive optical imaging with lanthanide lumiphores.

In principle, the millisecond emission lifetimes of lanthanide chelates should enable their ultrasensitive detection in biological systems by time-resolved optical microscopy. In practice, however, lanthanide imaging techniques have provided no better sensitivity than conventional fluorescence microscopy. Here, we identified three fundamental problems that have impeded lanthanide microscopy: low photon flux, inefficient excitation, and optics-derived background luminescence. We overcame these limitations with a new lanthanide imaging modality, transreflected illumination with luminescence resonance energy transfer (trLRET), which increases the time-integrated signal intensities of lanthanide lumiphores by 170-fold and the signal-to-background ratios by 75-fold. We demonstrate that trLRET provides at least an order-of-magnitude increase in detection sensitivity over that of conventional epifluorescence microscopy when used to visualize endogenous protein expression in zebrafish embryos. We also show that trLRET can be used to optically detect molecular interactions in vivo. trLRET promises to unlock the full potential of lanthanide lumiphores for ultrasensitive, autofluorescence-free biological imaging.

Determination of the conformational ensemble of the TAR RNA by X-ray scattering interferometry.

The conformational ensembles of structured RNA's are crucial for biological function, but they remain difficult to elucidate experimentally. We demonstrate with HIV-1 TAR RNA that X-ray scattering interferometry (XSI) can be used to determine RNA conformational ensembles. X-ray scattering interferometry (XSI) is based on site-specifically labeling RNA with pairs of heavy atom probes, and precisely measuring the distribution of inter-probe distances that arise from a heterogeneous mixture of RNA solution structures. We show that the XSI-based model of the TAR RNA ensemble closely resembles an independent model derived from NMR-RDC data. Further, we show how the TAR RNA ensemble changes shape at different salt concentrations. Finally, we demonstrate that a single hybrid model of the TAR RNA ensemble simultaneously fits both the XSI and NMR-RDC data set and show that XSI can be combined with NMR-RDC to further improve the quality of the determined ensemble. The results suggest that XSI-RNA will be a powerful approach for characterizing the solution conformational ensembles of RNAs and RNA-protein complexes under diverse solution conditions.

Directed Chemical Evolution with an Outsized Genetic Code.

The first demonstration that macromolecules could be evolved in a test tube was reported twenty-five years ago. That breakthrough meant that billions of years of chance discovery and refinement could be compressed into a few weeks, and provided a powerful tool that now dominates all aspects of protein engineering. A challenge has been to extend this scientific advance into synthetic chemical space: to enable the directed evolution of abiotic molecules. The problem has been tackled in many ways. These include expanding the natural genetic code to include unnatural amino acids, engineering polyketide and polypeptide synthases to produce novel products, and tagging combinatorial chemistry libraries with DNA. Importantly, there is still no small-molecule analog of directed protein evolution, i.e. a substantiated approach for optimizing complex (≥ 10^9 diversity) populations of synthetic small molecules over successive generations. We present a key advance towards this goal: a tool for genetically-programmed synthesis of small-molecule libraries from large chemical alphabets. The approach accommodates alphabets that are one to two orders of magnitude larger than any in Nature, and facilitates evolution within the chemical spaces they create. This is critical for small molecules, which are built up from numerous and highly varied chemical fragments. We report a proof-of-concept chemical evolution experiment utilizing an outsized genetic code, and demonstrate that fitness traits can be passed from an initial small-molecule population through to the great-grandchildren of that population. The results establish the practical feasibility of engineering synthetic small molecules through accelerated evolution.

The solution structural ensembles of RNA kink-turn motifs and their protein complexes.

With the growing number of crystal structures of RNA and RNA-protein complexes, a critical next step is understanding the dynamic solution behavior of these entities in terms of conformational ensembles and energy landscapes. To this end, we have used X-ray scattering interferometry (XSI) to probe the ubiquitous RNA kink-turn motif and its complexes with the canonical kink-turn binding protein L7Ae. XSI revealed that the folded kink-turn is best described as a restricted conformational ensemble. The ions present in solution alter the nature of this ensemble, and protein binding can perturb the kink-turn ensemble without collapsing it to a unique state. This study demonstrates how XSI can reveal structural and ensemble properties of RNAs and RNA-protein complexes and uncovers the behavior of an important RNA-protein motif. This type of information will be necessary to understand, predict and engineer the behavior and function of RNAs and their protein complexes.

From a structural average to the conformational ensemble of a DNA bulge.

Direct experimental measurements of conformational ensembles are critical for understanding macromolecular function, but traditional biophysical methods do not directly report the solution ensemble of a macromolecule. Small-angle X-ray scattering interferometry has the potential to overcome this limitation by providing the instantaneous distance distribution between pairs of gold-nanocrystal probes conjugated to a macromolecule in solution. Our X-ray interferometry experiments reveal an increasing bend angle of DNA duplexes with bulges of one, three, and five adenosine residues, consistent with previous FRET measurements, and further reveal an increasingly broad conformational ensemble with increasing bulge length. The distance distributions for the AAA bulge duplex (3A-DNA) with six different Au-Au pairs provide strong evidence against a simple elastic model in which fluctuations occur about a single conformational state. Instead, the measured distance distributions suggest a 3A-DNA ensemble with multiple conformational states predominantly across a region of conformational space with bend angles between 24 and 85 degrees and characteristic bend directions and helical twists and displacements. Additional X-ray interferometry experiments revealed perturbations to the ensemble from changes in ionic conditions and the bulge sequence, effects that can be understood in terms of electrostatic and stacking contributions to the ensemble and that demonstrate the sensitivity of X-ray interferometry. Combining X-ray interferometry ensemble data with molecular dynamics simulations gave atomic-level models of representative conformational states and of the molecular interactions that may shape the ensemble, and fluorescence measurements with 2-aminopurine-substituted 3A-DNA provided initial tests of these atomistic models. More generally, X-ray interferometry will provide powerful benchmarks for testing and developing computational methods.

Highly parallel translation of DNA sequences into small molecules.

A large body of in vitro evolution work establishes the utility of biopolymer libraries comprising 10(10) to 10(15) distinct molecules for the discovery of nanomolar-affinity ligands to proteins. Small-molecule libraries of comparable complexity will likely provide nanomolar-affinity small-molecule ligands. Unlike biopolymers, small molecules can offer the advantages of cell permeability, low immunogenicity, metabolic stability, rapid diffusion and inexpensive mass production. It is thought that such desirable in vivo behavior is correlated with the physical properties of small molecules, specifically a limited number of hydrogen bond donors and acceptors, a defined range of hydrophobicity, and most importantly, molecular weights less than 500 Daltons. Creating a collection of 10(10) to 10(15) small molecules that meet these criteria requires the use of hundreds to thousands of diversity elements per step in a combinatorial synthesis of three to five steps. With this goal in mind, we have reported a set of mesofluidic devices that enable DNA-programmed combinatorial chemistry in a highly parallel 384-well plate format. Here, we demonstrate that these devices can translate DNA genes encoding 384 diversity elements per coding position into corresponding small-molecule gene products. This robust and efficient procedure yields small molecule-DNA conjugates suitable for in vitro evolution experiments.

Mesofluidic devices for DNA-programmed combinatorial chemistry.

Hybrid combinatorial chemistry strategies that use DNA as an information-carrying medium are proving to be powerful tools for molecular discovery. In order to extend these efforts, we present a highly parallel format for DNA-programmed chemical library synthesis. The new format uses a standard microwell plate footprint and is compatible with commercially available automation technology. It can accommodate a wide variety of combinatorial synthetic schemes with up to 384 different building blocks per chemical step. We demonstrate that fluidic routing of DNA populations in the highly parallel format occurs with excellent specificity, and that chemistry on DNA arrayed into 384 well plates proceeds robustly, two requirements for the high-fidelity translation and efficient in vitro evolution of small molecules.

Structural and kinetic mapping of side-chain exposure onto the protein energy landscape.

Identification and characterization of structural fluctuations that occur under native conditions is crucial for understanding protein folding and function, but such fluctuations are often rare and transient, making them difficult to study. Native-state hydrogen exchange (NSHX) has been a powerful tool for identifying such rarely populated conformations, but it generally reveals no information about the placement of these species along the folding reaction coordinate or the barriers separating them from the folded state and provides little insight into side-chain packing. To complement such studies, we have performed native-state alkyl-proton exchange, a method analogous to NSHX that monitors cysteine modification rather than backbone amide exchange, to examine the folding landscape of Escherichia coli ribonuclease H, a protein well characterized by hydrogen exchange. We have chosen experimental conditions such that the rate-limiting barrier acts as a kinetic partition: residues that become exposed only upon crossing the unfolding barrier are modified in the EX1 regime (alkylation rates report on the rate of unfolding), while those exposed on the native side of the barrier are modified predominantly in the EX2 regime (alkylation rates report on equilibrium populations). This kinetic partitioning allows for identification and placement of partially unfolded forms along the reaction coordinate. Using this approach we detect previously unidentified, rarely populated conformations residing on the native side of the barrier and identify side chains that are modified only upon crossing the unfolding barrier. Thus, in a single experiment under native conditions, both sides of the rate-limiting barrier are investigated.

Design of protein-ligand binding based on the molecular-mechanics energy model.

While the molecular-mechanics field has standardized on a few potential energy functions, computational protein design efforts are based on potentials that are unique to individual laboratories. Here we show that a standard molecular-mechanics potential energy function without any modifications can be used to engineer protein-ligand binding. A molecular-mechanics potential is used to reconstruct the coordinates of various binding sites with an average root-mean-square error of 0.61 A and to reproduce known ligand-induced side-chain conformational shifts. Within a series of 34 mutants, the calculation can always distinguish between weak (K(d)>1 mM) and tight (K(d)<10 microM) binding sequences. Starting from partial coordinates of the ribose-binding protein lacking the ligand and the 10 primary contact residues, the molecular-mechanics potential is used to redesign a ribose-binding site. Out of a search space of 2 x 10(12) sequences, the calculation selects a point mutant of the native protein as the top solution (experimental K(d)=17 microM) and the native protein as the second best solution (experimental K(d)=210 nM). The quality of the predictions depends on the accuracy of the generalized Born electrostatics model, treatment of protonation equilibria, high-resolution rotamer sampling, a final local energy minimization step, and explicit modeling of the bound, unbound, and unfolded states. The application of unmodified molecular-mechanics potentials to protein design links two fields in a mutually beneficial way. Design provides a new avenue for testing molecular-mechanics energy functions, and future improvements in these energy functions will presumably lead to more accurate design results.

Synthetic ligands discovered by in vitro selection.

The recognition and catalytic properties of biopolymers derive from an elegant evolutionary mechanism, whereby the genetic material encoding molecules with superior functional attributes survives a selective pressure and is propagated to subsequent generations. This process is routinely mimicked in vitro to generate nucleic-acid or peptide ligands and catalysts. Recent advances in DNA-programmed organic synthesis have raised the possibility that evolutionary strategies could also be used for small-molecule discovery, but the idea remains unproven. Here, using DNA-programmed combinatorial chemistry, a collection of 100 million distinct compounds is synthesized and subjected to selection for binding to the N-terminal SH3 domain of the proto-oncogene Crk. Over six generations, the molecular population converges to a small number of novel SH3 domain ligands. Remarkably, the hits bind with affinities similar to those of peptide SH3 ligands isolated from phage libraries of comparable complexity. The evolutionary approach has the potential to drastically simplify and accelerate small-molecule discovery.

Chemical evolution as a tool for molecular discovery.

In modern academic and industrial laboratories, evolutionary strategies are used routinely to identify biopolymers with novel activities. Large libraries of nucleic acids (approximately 10(15)) or peptides and proteins (approximately 10(13)) can be subjected to multiple rounds of selective pressure, amplification, and diversification, yielding individual sequences with desirable properties. Although the evolutionary approach is a powerful search tool, the chemical nature of biopolymers is not suited for all purposes. Application of evolutionary strategies to libraries of arbitrary chemical composition would overcome this problem, and radically change traditional small-molecule discovery. The chemical make-up of in vitro evolution libraries has necessarily been limited, because library synthesis relies on enzymes. A great deal of current research focuses on expanding the chemical repertoire of in vitro evolution by (a) broadening enzyme substrate specificities to include unnatural building blocks, or (b) developing methods to translate DNA sequences into multistep organic syntheses. We discuss the strengths and weaknesses of the approaches, review the successes, and consider the future of chemical evolution as a tool.

Potential energy functions for protein design.

Different potential energy functions have predominated in protein dynamics simulations, protein design calculations, and protein structure prediction. Clearly, the same physics applies in all three cases. The differences in potential energy functions reflect differences in how the calculations are performed. With improvements in computer power and algorithms, the same potential energy function should be applicable to all three problems. In this review, we examine energy functions currently used for protein design, and look to the molecular mechanics field for advances that could be used in the next generation of design algorithms. In particular, we focus on improved models of the hydrophobic effect, polarization and hydrogen bonding.

Misincorporation proton-alkyl exchange (MPAX): engineering cysteine probes into proteins.

This unit describes a rapid and efficient method to screen a polypeptide for amino acid residues that contribute to protein-protein interaction interfaces. Cysteine residues are introduced as positional probes in a protein at random by co-expression in bacteria with specific cysteine misincorporator tRNAs. The protein is then purified as an ensemble of polypeptides containing cysteine at low frequency, at different positions in each molecule. The ability of the native protein structure to protect different cysteine residues from chemical modification by iodoacetamide is determined to obtain a protein surface map that reveals candidate surface residues that are likely to be important for protein-protein interaction. Cysteine mutants with altered ligand binding can also be selected simultaneously by affinity chromatography.

Structural test of the parameterized-backbone method for protein design.

Designing new protein folds requires a method for simultaneously optimizing the conformation of the backbone and the side-chains. One approach to this problem is the use of a parameterized backbone, which allows the systematic exploration of families of structures. We report the crystal structure of RH3, a right-handed, three-helix coiled coil that was designed using a parameterized backbone and detailed modeling of core packing. This crystal structure was determined using another rationally designed feature, a metal-binding site that permitted experimental phasing of the X-ray data. RH3 adopted the intended fold, which has not been observed previously in biological proteins. Unanticipated structural asymmetry in the trimer was a principal source of variation within the RH3 structure. The sequence of RH3 differs from that of a previously characterized right-handed tetramer, RH4, at only one position in each 11 amino acid sequence repeat. This close similarity indicates that the design method is sensitive to the core packing interactions that specify the protein structure. Comparison of the structures of RH3 and RH4 indicates that both steric overlap and cavity formation provide strong driving forces for oligomer specificity.

DNA display I. Sequence-encoded routing of DNA populations.

Recently reported technologies for DNA-directed organic synthesis and for DNA computing rely on routing DNA populations through complex networks. The reduction of these ideas to practice has been limited by a lack of practical experimental tools. Here we describe a modular design for DNA routing genes, and routing machinery made from oligonucleotides and commercially available chromatography resins. The routing machinery partitions nanomole quantities of DNA into physically distinct subpools based on sequence. Partitioning steps can be iterated indefinitely, with worst-case yields of 85% per step. These techniques facilitate DNA-programmed chemical synthesis, and thus enable a materials biology that could revolutionize drug discovery.